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We study an ensemble of random walkers carrying internal noisy phase oscillators which are 
synchronized among the walkers by local interactions. Due to individual mobility, the interaction 
partners of every walker change randomly, hereby introducing an additional, independent source of 
fluctuations, thus constituting the intrinsic nonequilibrium nature of the temporal dynamics. We 
employ this paradigmatic model system to discuss how the emergence of order is affected by motion 
of individual entities. In particular, we consider both, normal diffusive motion and superdiffusion. 
A non-Hamiltonian field theory including multiplicative noise terms is derived which describes the 
nonequilibrium dynamics at the macroscale. This theory reveals a defect-mediated transition from 
incoherence to quasi long-range order for normal diffusion of oscillators in two dimensions, implying 
a power-law dependence of all synchronization properties on system size. In contrast, superdif- 
fusive transport suppresses the emergence of topological defects, thereby inducing a continuous 
synchronization transition to long-range order in two dimensions. These results are consistent with 
particle-based simulations. 
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I. INTRODUCTION 

In a large variety of chemical and biological systems, 
moving entities synchronize an internal degree of free¬ 
dom: examples include mobile catalytic beads exhibiting 
the oscillatory Belousov-Zhabotinsky (BZ) reaction [T] 
as well as the BZ reaction in vortical flow [2], synthetic 
genetic oscillators [5] , motile cells synchronizing intracel¬ 
lular oscillations during embryonic somitogenesis gHS], 
quorum sensing of signaling amoebae [7] and the coor¬ 
dinated motion of myxobacteria regulated by the oscil¬ 
latory Frz-signaling system [HHin]. Furthermore, it has 
been speculated that the flagellar beating of active mi¬ 
croswimmers synchronizes by hydrodynamic interaction 
thus affecting the rheology of active suspensions mm- 

Motivated by the above-mentioned nonequibrium sys¬ 
tems, we study synchronization in an ensemble of ei¬ 
ther diffusively or superdiffusively moving random walk¬ 
ers carrying internal noisy oscillators which synchronize 
by local interactions. This paradigmatic model describes 
a nonequilibrium system as evidenced by an analogy 
to Monte-Carlo algorithms m- random motion corre¬ 
sponds to Kawasaki dynamics m at infinite tempera¬ 
ture, whereas the oscillator dynamics in the presence of 
noise is similar to Glauber updates m at a finite tem¬ 
perature. Ergo, the system is coupled to two heat baths 
simultaneously and, as a consequence, does never reach 
thermodynamic equilibrium. 

We are primarily addressing the question how the 
emergence of order in spatially extended systems, which 
relies on the suppression of fluctuations by sufficiently 
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fast information transfer through the system, is influ¬ 
enced by the motion type of individual entities. This 
question is of interest beyond synchronization in nonequi¬ 
librium statistical mechanics, ranging from active mat¬ 
ter [Mm and collective motion [B [2D] to stochastic 
reaction-diffusion systems mm and epidemic spread¬ 
ing [D51 - [D5] . Ultimately, the interplay of local interaction 
and transport resulting in the emergence of order is a key 
element in all the above-named applications. 

The study of phase transitions and spontaneous sym¬ 
metry breaking in isothermal systems at thermodynamic 
equilibrium has already led to fundamental theoretical 
insights regarding the emergence of order such as the uni¬ 
versality of critical behavior |2D| or the Mermin-Wagner 
theorem 1211: which excludes the emergence of long-range 
order in low-dimensional equilibrium systems with a con¬ 
tinuous symmetry at finite temperature. This work fo¬ 
cusses on the crucial influence of the motility of indi¬ 
vidual entities on the emergence of macroscopic order in 
nonequilibrium systems, which is less well understood. 

The temporal dynamics of many oscillators models, 
such as phase oscillators which constitute one of the cor¬ 
ner stones of synchronization theory |28j , is invariant un¬ 
der global phase shifts. This leads to the emergence of 
spin-waves in one dimension IBIDD] and topological de¬ 
fects [3T] in two dimensions. These facts, that relate 
synchronization theory and statistical mechanics, were 
often overlooked in studies of motile oscillators [321438] . 
Here, we discuss how the motion of these topological ex¬ 
citations and the corresponding annihilation dynamics is 
affected by the motion of the oscillators. 

Whereas various previous studies focused on determin¬ 
istic oscillators [32ll36| . we examine fluctuation-induced 
properties of synchronized states such as correlation func¬ 
tions and study synchronization from a statistical me- 
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chanics perspective. In particular, we systematically de¬ 
rive a field theory from the oscillator model and demon¬ 
strate that oscillators moving diffusively exhibit a defect- 
mediated transition from incoherence to quasi long-range 
order (QLRO) analogous to the Berezinskii-Kosterlitz- 
Thouless (BKT) transition [3^ 00] in two dimensions. 
We show further that a continuous nonequilibrium tran¬ 
sition to long-range order (LRO) does occur - forbid¬ 
den in equilibrium systems according to the Mermin- 
Wagner theorem [27] - if oscillators move superdiffu- 
sively, thereby accelerating the annihilation of defects. 


II. MODEL 


We consider an ensemble of N random walkers in a 
periodic space, each of them carrying an internal clock. 
Individual clocks, which are assumed to possess identical 
mean natural frequencies wq, are modeled as noisy phase 
oscillators that synchronize to the clocks of neighboring 
walkers. The state of the j-th walker is thus characterized 
by its position rj{t) and the phase 9j{t) or the relative 
phase Xj{t) = ^j{i) ~ ^ot- The dynamics of rj(t) and 
Xjit) obeys the Langevin equations 

rj = |j(t), (la) 

= - rj\)sm{xk - Xj) + \/^^Vjit)-{^'^) 


This model system can be understood as a nonequilib¬ 
rium, off-lattice version of the XY model where individ¬ 
ual spins are motile. 

The random motion of particles (or spins) is de¬ 
scribed by Eq. (la), where ^j{t) denotes Levy noise [HI 
02] ■ Hence, a particle jumps according to rj{t + T) = 
rj{t) + dt' in a time interval t > 0. The 

random displacements follow a Levy a-stable distribu¬ 
tion 03] defined by its propagator Ga{i^,T) in Fourier 
domain: (^^(q,r) = . The motion is character¬ 

ized by a generalized diffusion coefficient 7 determining 
the width of the displacement distribution, and an in¬ 
dex a describing the tails of the step length distribution. 
For a = 2, the displacement distribution is Gaussian and, 
consequently, particles undergo Brownian motion with a 
diffusion coefficient 7 . In contrast, the step length distri¬ 
bution is heavy tailed for aG ( 0 , 2 ) implying large jumps. 
Hence, walkers perform Levy flights mM- 

Eq. (lb I describes the oscillator dynamics. The first 
term accounts for the local interaction among oscillators 
analogous to the paradigmatic Kuramoto model [HI 051 - 
03 . The interaction strength /? of two oscillators is a 
function of their relative distance. We assume that the in¬ 
teraction is strictly short-ranged, i.e. it vanishes beyond a 
certain interaction range representing the sensing radius 
of an individual oscillator that we set to one: /3{x) = 0 
for |x| > 1. The interaction is rescaled by the number 
of neighbors Nj = ~ I’il) reflecting that the 

adaptivity of an oscillator is independent of the number 


of neighbors. The second term in Eq. is an additive, 
Gaussian white noise |48| with intensity D^. 


III. PHENOMENOLOGY 


We performed numerical simulations particularly 
studying the influence of oscillator motion on synchro¬ 
nization in two dimensions. Here, the density of particles 
is set to po = 1 ; which is below the percolation thresh- 
old Po.perc. ~ 1.44 m- Therefore, the communication 
of oscillators over large scales is dominated by particle 
transport. 

First, we consider the order parameter $ = 








as a function of the noise strength 


(left column in Fig. [^. A transition from incoherence 
(4) ~ 0 ) to synchronization ($ > 0 ) is observed below 
a critical noise strength The dependence of 

on the diffusivity 7 indicates that synchronization is en¬ 
hanced by increasing 7 . For 7 —>■ 00 , the order parameter 
tends to the order parameter for globally coupled oscil¬ 
lators 03]. We can conclude that motion enhances in¬ 
formation transfer in finite systems, thus promoting the 
emergence of order. 

The transport properties crucially influence the syn¬ 
chronization properties which becomes evident from 
varying the particle number N and system size L such 
that the average density po = N /is constant (finite- 
size scaling). For Brownian motion, the order parameter 
decreases according to a power law with the particle num- 
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FIG. 1. Order parameter 4' (averaged over time) as a function 
of the noise (left) and the particle number N (right) for 
oscillators moving diffusively (upper row) and by Levy flights 
(lower row). The blue dashed line represents the solution of 
the Kuramoto model (critical point: D® = 1/2) for global 
coupling 09]. In accordance with the theory (cf. section IV I, 
the lines in the upper right panel correspond to power law 
fits, while the horizontal lines in the lower right panel indi¬ 
cate the saturation of the order parameter with system size. 
Parameters: po = 1, Pix) = 0(1 — x), L = 300 (left column), 
7 = 0.1 (right column), numerical time step r = 0 . 01 . 
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FIG. 2. Correlation function (averaged over time) of s(xi) = 
(cos Xi, sinxO as a function of the distance for diffusing oscil¬ 
lators. For high noise (disordered phase), correlations decay 
exponentially (left panel). In contrast, correlations decrease 
algebraically for low noise (right panel) confirming the BKT 
nature of the synchronization transition and the emergence 
of QLRO. Lines represent fits using the associated functional 
form. Parameters: a = 2, 7 = 0.1, po = 1, L = 300, 
I3{x) = 0(1 — x), numerical time step r = 0.01. 


her (top right in Fig.[^. The corresponding exponent is a 
function of the model parameters, particularly the noise 
strength D^. Thus, the system is globally disordered in 
the thermodynamic limit. Similar information is con¬ 
veyed by the correlation function (Fig. [^. We assigned 
to each phase variable Xi ^ spin s(xi) = (cosxi, sinxi) 
and calculated the correlations of s as a function of 
the distance. Correlations decay exponentially for high 
noise but decrease according to a power law in the low 
noise regime as QLRO emerges through a BKT transi¬ 
tion [Ml H(T| . 

Levy flights change the qualitative picture signifi¬ 
cantly: the order parameter $ tends towards a nonzero 
constant value for large system sizes (bottom right in 
Fig. 0 suggesting the emergence of a robust synchronized 
state. Further, this result implies, notably, the emergence 
of long-range order (LRO) in two dimensions. We find 
that superdiffusion, in contrast to what we observe for 
diffusive transport, qualitatively affects the large-scale 
properties of the ordered states and, consequently, the 
nature of the associated phase transition. In section [TV] 
we confirm these findings by a field theoretic analysis. 
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FIG. 3. Relaxational dynamics from random initial condi¬ 
tions: vortices (V) and anti-vortices (AV) are schematically 
shown by vectors s(xj) = (cosxj, sinxi) in panels (a,b) and 
(c,d), respectively. A snapshot of a vortex-anti-vortex pair is 
depicted in (e), where phases Xj nre color-coded. 



FIG. 4. Relaxational dynamics from random initial condi¬ 
tions: defect density n{t) (averaged over multiple realizations) 
for Brownian motion (top) and Levy flights (bottom) as a 
function of time. On the right, the defect density is rescaled 
by t/\nt. The black line in (a) also indicates this scaling. 
Superdiffusive transport (d) qualitatively enhances the relax¬ 
ational dynamics as compared to the case of normal diffusion. 
Parameters: po = L = 0, L — 300 (a,b), L = 400 (c,d), 
P{x) = 0(1 — x), numerical time step r = 0 . 01 . 


The emergence of long-range order in the case of su¬ 
perdiffusive transport can intuitively be understood by 
studying in which way macroscopic fluctuations decay. 
We quantify the temporal relaxation to the stationary 
state by preparing the system in the disordered state 
and, thereafter, observing how the system relaxes for 
zero noise = 0). After an initial transient, the sys¬ 
tem is locally ordered but globally disordered due to the 
emergence of free topological defects (vortices and anti¬ 
vortices) which eventually annihilate pairwise (see Fig. 
for a visualization). In general, defects annihilate faster if 
the diffusivity of oscillators is high (cf. Fig. [^. However, 
the asymptotic behavior of the time dependent defect 
density n(t) does not depend on the diffusion coefficient 7 
for normal diffusion (black line in Fig. 1^): n{t) ~ Int/t 
similar to the classical XY model [51 j7 On the other 
hand, the annihilation process is qualitatively acceler¬ 
ated towards a power-law like decay of the defect density 
n{t) ^ with (5 > 1 in the long-time limit, if particles 
perform Levy flights indicating that superdiffusion of in¬ 
dividual entities induces a novel regime of the relaxation 
process. 


IV. FIELD THEORY 

We corroborate the simulation results by deriving and 
analyzing field equations for the relevant order parame¬ 
ters. The instantaneous particle density 

pir,t)=j:US{r-r,it)) ( 2 ) 

is a central quantity of interest, since the presence of den¬ 
sity fluctuations distinguishes the nonequilibrium tempo¬ 
ral dynamics of motile oscillators from lattice systems in 
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equilibrium - the dynamic network of motile particles in¬ 
duces a fluctuating number of neighbors thus providing a 
source of nonthermal noise. In the stationary state, the 
mean density is constant since there are no forces act¬ 
ing among the particles thus excluding the emergence of 
density instabilities. Fluctuations around the mean are 
determined by the propagator (r, t) which also reflects 
the transition probability for a spatial displacement of a 
particle within a certain time interval: 


/d-^r/JdrDIrl 


^ /\AdJd'^r'f}{\r'\)\ stems from the inter¬ 


action among the oscillators; the term in brackets indi¬ 
cates the emergence of local order via a pitchfork bifurca¬ 
tion for = 1/2, where e = 1/2. Fluctuations on 

the macroscale are reflected by the Gaussian field ri{r,t) 
in Eq. Q, which is projected by the matrix 


n = 




^Jl-\S/p\'^MW. 


{p{r,t)) = po (3a) 

(p(r, t)p{r', t')) - pI = pq [G„(r t') - . (3b) 

Correlations of the density field are proportional to the 
mean density and, therefore, the variance of the rescaled 
density Sp = [p/po — 1] decreases proportional to Pq^ in 
accordance with the central limit theorem |52j . 

Along with the density describing the spatial distribu¬ 
tion of particles, we introduce an order parameter field 


S(r,t) 


|S| 


/cos'(/'\ /cos Xj{t) 

^ l^sinxi(t) 


(5(r 


(4) 


reflecting the local degree of synchronization. This field 
variable is related to the global order parameter $ that 
was used to characterize particle based simulations via 



^y'd"rdV(S(r,t).S(r',t)). 


Hence, an increase of the global order parameter - as 
observed in numerical simulations - is related to a quali¬ 
tative change of the correlation function of S: correlations 
become long-ranged at the critical point. 

The definition of the order parameter, Eq. Q, con¬ 
stitutes a change of variables from {rj(t),Xj(t)} to the 
field S(r, t). Therefore, the dynamics is obtained by ap¬ 
plication of stochastic calculus [3S1 [S3] . It turns out that 
the dynamics is coupled to the density and the nematic 
tensor Q = Ylf=i-^[Xj{'t)] '^(r “ where 

_ /^cos(2(^) sin{2ip)\ 

^ “ ^sin {2p) - cos {2p)) ■ 

Here, we use the closure relation |S|^M[?/>]/p to ob¬ 
tain a self-consistent field equation for S(r, t): 


5tS(r,t) « 7 A“/ 2 S-kMS + 




+ \jD^p{Y,t)/2 n • T7(r, t). 


S 

( 6 ) 


This field theory reflects the microscopic dynam¬ 
ics: the (fractional) Laplacian multiplied by the diffu- 
sivity 7 describes the transport of particles [411 HI] 
giving rise to (non-)local coupling; the coefficient /3 = 


We emphasize (i) the multiplicative noise comprising the 
square root of the density and (ii) the non-Hamiltonian 
dynamics. Consequently, the probability to find a cer¬ 
tain field configuration is generally not determined by a 
Boltzmann distribution, in turn reflecting the absence of 
temperature and the nonequilibrium nature of the model. 
In this regard, the dynamics of the system close to the 
transition point may potentially contain novel physics 
including a set of non-classical critical exponents in the 
low density regime. These questions which go beyond 
the scope of this work may be addressed in future stud¬ 
ies. Here, we focus on the properties of the synchronized 
phase as well as the emergent dynamics of topological 
defects. 

In the high density regime, fluctuations of the density 
can approximately be neglected according to Eq. ([^ en¬ 
abling us to link the nonequilibrium system featuring a 
dynamic network of short-ranged interactions to the field 
theory of the static XY model [51H57]. This nontriv¬ 
ial mapping, which connects two fundamentally different 
physical regimes, explains the phenomenology observed 
in simulations in several aspects. We deduce from Eq. ([^ 
that the diffusivity 7 rescales the coupling coefficient (3 
in the case of Brownian diffusion (a = 2). The struc¬ 
ture of Eq. however, does not change for 7 —>• 0 in 
that case. Therefore, the qualitative phenomenology of 
the XY model is observed if oscillators perform Brow¬ 
nian motion including a BKT transition towards QLRO 
whereas the equation is structurally different for superdif- 
fusive motion. 

Close to the critical point (in the disordered phase), 
the local order parameter is small, |S//o| ^ 1, and there¬ 
fore n « 2 • 1 to lowest order. In this limit, the field 
equation ([^ reduces to the Ginzburg-Landau model with 
additive white noise enabling us to determine the quali¬ 
tative dependence of the transition point on model 
parameters [331133] : 


Xc Xc 


1 -^ 




dy 


PoJ (27r)‘'7|q|“-|-,5|q|' 


0{e^) 


■ (7) 


The integral in Eq. Q runs over all wavevectors 
with |q| < A, where A is a cutoff parameter inversely 
proportional to the coarse-graining scale. Since the inte¬ 
grand is positive, the critical noise is shifted towards 

smaller noise values with respect to D^l as observed in 
simulations (cf. Eig.[^. Only if the mobility or the parti¬ 
cle density is high, the noise-induced shift of vanishes 
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and the order-disorder transitions on the local and global 
level occur simultaneously: ~ D^xl- 

We focus now on the properties of locally ordered 
states, where S(r,t) « S'o(cos'!/'(i',t),sin'0(r,t)). In 
this regime, all quantities are determined by the 
phase field such as the correlation func¬ 

tion (S(r,<)-S(r',t)) « (cos <) — ?/>(r', t)]). The 
dynamics of ■(/’(r, t) follows from Eq. §: 


dt^ 


Q f^d^, sin[^(r',f) - V>(r,t)] 

J |r' - 




The noise strength is denoted by >f = D^pq{1 -|- 
Sq/pq)/Sq, C is a constant [60] and <;(r,t) is a 
Gaussian white field m- We assume further 
that tp(r,t) varies slowly such that we can approximate 
sin['!/;(r', t) — t)] « t) — t)] to obtain 


dt'ilj{r,t) ^ + jSA-ip + ( 8 ) 


The solution of this linear equation is straightforward via 
Fourier transformation. One finds that the stationary 
phase field is Gaussian allowing us to calculate 

the order parameter correlations by the identity m 


(S(r,t)-S(r',t))«S'^exp 




Accordingly, the order parameter correlation function is 
determined by the mean squared phase difference which 
follows from the solution of Eq. (|^ : 


(bP{r,t) 


= X 


(27ry 7|q|“-h/3|q|^ 


(9) 


The mean squared phase difference at two points r and r' 
behaves for large separations as |r — i.e. it decays 

for a < d indicating a long-ranged ordered state due to 
superdiffusive motion of particles. Note that the large 
scale dynamics essentially depends on the dimension of 
space. In contrast, microscopic details do not change the 
large scale properties qualitatively, thus indicating a form 
of universality. 

Furthermore, the deterministic equation for ^/)(r,t) 
supports vortex solutions V^(r,t) = ± arg(r) -|- tpo in 
two dimensions as observed in numerical simulations 
(cf. Fig.^lU). We exploit Eq. ^ to determine the relax- 
ational dynamics of the defect density from a qualitative 
scaling argument. By successive rescaling of space r — 6 r 
and time t —>■ b°‘t with a positive scale parameter &, the 
diffusive coupling {f3 —)■ as well as the noise term 

(>!:—>■ b°‘~'^ k) are rendered irrelevant for a < 2 whereas 7 
remains unchanged. Thus, the phase equation becomes 
scale invariant in the long time limit. Hence, the mean 
defect number in a cube of volume at time t is asymp¬ 
totically equal to the defect number in the domain of size 
{bL^ at time b°‘t (cf. Fig. ^ for an illustration). Conse¬ 
quently, the defect density n{t) is asymptotically a homo¬ 
geneous function decaying as n(t) oc where the ex¬ 

ponent explicitly depends on the motion type via a, but is 
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FIG. 5. Illustration of the self-similarity argument determin¬ 
ing the decay of the density of topological defects. The dy¬ 
namics of the phase field cf. Eq. (j^, is asymptoti¬ 

cally scale invariant under the rescaling r — >■ 6r and t —>■ 6“t. 
Therefore, the mean number of defects M in a cube of size 
at time t (left panel) and in a cube (bL)^ at time b°‘t (right 
panel) is approximately equal, as illustrated above. 


independent of microscopic details. Accordingly, the re¬ 
laxation is accelerated in the case of superdiffusive trans¬ 
port as observed in numerical simulations (cf. Fig. |^. 


V. SUMMARY & OUTLOOK 

We studied an ensemble of random walkers carrying 
internal clocks which are synchronized by local interac¬ 
tions. The phenomenology of this paradigmatic nonequi¬ 
librium model, as observed in numerical simulations, is 
remarkably rich. Numerical observations were corrobo¬ 
rated by a non-Hamiltonian field theory including multi¬ 
plicative noise terms which was analytically derived from 
the microscopic particle dynamics. The field theory re¬ 
veals that the qualitative phenomenology of the system 
is not determined by mixing rates or, equivalently, the 
diffusion coefficients but rather by the motion type of 
particles. 

For Brownian diffusion of individual oscillators, we re¬ 
port on a Berezinskii-Kosterlitz-Thouless transition from 
incoherence to quasi long-range order accompanied by 
the emergence of topological defects implying power-law 
like system size dependencies of all synchronization prop¬ 
erties. We argued that the nonequilibrium system with 
diffusive particle transport exhibits the same qualitative 
behavior as an equilibrium system where particles are 
held fixed on a lattice. This implies that long-range or¬ 
der cannot emerge for dimensions smaller or equal to two 
in the thermodynamic limit and, thus, extends the pre¬ 
dictions of the Mermin-Wagner theorem m- valid for 
equilibrium systems - to a nonequilibrium system with 
diffusive particle transport. 

Furthermore, we showed that superdiffusive transport 
may give rise to long-range order in two and even in one 
dimension depending on the characteristics of the mo¬ 
tion of individual entities. The emergence of long-range 
order in two dimensions relies on the accelerated anni¬ 
hilation dynamics of topological defects as confirmed by 
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numerical simulations in two dimensions. 

In short, we characterized in this work the condi¬ 
tions required for the emergence of long-range order in 
nonequilibrium systems. We believe that the stochastic 
dynamics of moving oscillators provides an interesting 
playground to study the interplay of local interaction and 
particle transport in the context of nonequilibrium statis¬ 
tical mechanics from a theoretical point of view, involv¬ 
ing the mixture of normal fluctuations (oscillator dynam¬ 
ics) and anomalous fluctuations (superdiffusive motion) 
as well as the emergence of order and dynamics of topo¬ 
logical defects in this context. Besides, the results may 
serve as a reference point for the description of specific 


experiments such as those mentioned in the introduction, 
where modeling may require a higher level of complexity, 
e.g. the addition of time delay |55] . 
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